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Abstract 

This work is devoted to prove new entropy estimates for a general isothermal 
model of capillary fluids derived by J.E Dunn and J.Serrin (1985) (see 53j)j which 
can be used as a phase transition model. More precisely we will derive new estimates 
for the density and we will give a new structure for the Korteweg system which 
allow us to obtain the existence of global weak solution. The key of the proof comes 
from the introduction of a new effective velocity. The proof is widely inspired from 
the works of A. Mellet and A. Vasseur (see [Ml)- In a second part, we shall give a 
Prody-Serrin blow-up criterion for this system which widely improves the results of 
|17] and the known results on compressible systems. 

1 Introduction 

We are concerned with compressible fluids endowed with internal capillarity. The model 
we consider originates from the XlXth century work by Van der Waals and Korteweg 
[4HI29j and was actually derived in its modern form in the 1980s using the second gradient 
theory, see for instance \14:\ l28l HQ] . The first investigations begin with the Young-Laplace 
theory which claims that the phases are separated by a hypersurface and that the jump 
in the pressure across the hypersurface is proportional to the curvature of the hypersur- 
face. The main difficulty consists in describing the location and the movement of the 
interfaces. 

Another major problem is to understand whether the interface behaves as a discontinu- 
ity in the state space (sharp interface) or whether the phase boundary corresponds to 
a more regular transition (diffuse interface, DI). The diffuse interface models have the 
advantage to consider only one set of equations in a single spatial domain (the density 
takes into account the different phases) which considerably simplifies the mathematical 
and numerical study (indeed in the case of sharp interfaces, we have to treat a problem 
with free boundary). 

Another approach corresponds to determine equilibrium solutions which classically con- 
sists in the minimization of the free energy functional. Unfortunately this minimization 
problem has an infinity of solutions, and many of them are physically wrong (some de- 
tails are given later). In order to overcome this difficulty, Van der Waals in the XlX-th 
century was the first to add a term of capillarity to select the physically correct solutions, 
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modulo the introduction of a diffuse interface. This theory is widely accepted as a ther- 
modynamically consistent model for equilibria. Alternatively, another way to penalize 
the high density variations consists in applying a zero order but non-local operator to 
the density gradient [36], [37], |38j . We refer for a mathematical analysis on this system 
to [121 US [22l [23]. Korteweg-type models are based on an extended version of nonequi- 
librium thermodynamics, which assumes that the energy of the fluid not only depends 
on standard variables but also on the gradient of the density. 

Let us now consider a fluid of density p > 0, velocity fleld u G Q.^ (both are deflned on a 
subset Vl with Vl = $7^ or the torus T^), we are now interested in the following compress- 
ible capillary fluid model, which can be derived from a Cahn-Hilliard like free energy (see 
the pioneering work by J.- E. Dunn and J. Serrin in [H] and also in [21 121 [TU| I27|). The 
conservation of mass and of momentum write: 



il.l) 



( d 

—p + div{pu) = 0, 
d 

^ dt^^^^ ~^ div(/9n u) — div(^p Vu) — div(ap Vu*) -|- \/{ap'^) = divK, 

where the Korteweg tensor read as following: 

divK = V {pK{p)Ap + ^{k{p) + pK {p))\V pI"^) - div {k{p)Vp0Vp). (1.2) 

K is the coefficient of capillarity and is a regular function of the form K.{p) = Kp" with 
a £ n. The term divi^ allows to describe the variation of density at the interfaces 
between two phases, generally a mixture liquid-vapor. P = ap^ with 7 > 1 is a general 
barotropic pressure term, p and a are the two Lame viscosity coefficients and satisfying: 

p > and p > a > 0. 

In particular, it allows to write the diffusion tensor on the form (p — a)div(pVn) -|- 
adiv{pDu) with Du = Vn-|-* Vn the strain tensor which implies some energy inequality. 
More precisely when we multiply the momentum equation by n, we have: 

{p{t,x)\u{t,x)\^ + ::^p'{t,x) + K\v^Wt,x))dx 

+ [ [ {{p - a)p{t,x)\Vu\'^{t,x) + ap{t,x)\Du\'^{t,x) (1.3) 
Jo Jn 

1 



< C{ / (po(x)|^;o(x)P + -pI{x) + K|VVp^|2)dx). 

./o 7 — i 



n 7 

We now want to recall the existing results on the existence of global weak solutions on 
classical compressible Navier Stokes equations and on Korteweg system. 



Weak solutions 

We wish to prove existence and uniqueness results for [NHV) in functions spaces very 
close to energy spaces. In the non isothermal non capillary case and P{p) = ap^ , with 
a > and 7 > 1, P-L. Lions in [32] proved the global existence of variational solutions 
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(p, u, 9) to (NHV) with k = for 7 > f if iV > 4, 7 > ^ if iV = 2, 3 and initial data 
{po,mo) such that: 

I 1 2 

n(po)-n(p), &eLn^^^), and poOoeL^n^). 
Po 

These solutions are weak solutions in the classical sense for the equation of mass conser- 
vation and for the equation of the momentum. Notice that the main difficulty for proving 
Lions' theorem consists in exhibiting strong compactness properties of the density p in 
-^foc spaces required to pass to the limit in the pressure term P{p) = ap^ . 
Let us mention that Feireisl in ^15j generalized the result to 7 > ^ i'^ establishing that we 
can obtain renormalized solution without imposing that p G Lf^,^-, for this he introduces 
the concept of oscillation defect measure evaluating the lost of compactness. 
We can finally cite the result from Bresch-Desjardins in [5], [6] where they show the exis- 
tence of global weak solution for [NHV) with k = and with a cold pressure. In fact 
they are working with specific viscosity coefficients which verify an algebraic relation. It 
allows them to get good estimate on the density by using new energy inequality and to 
treat by compactness all the delicate terms as the pressure. In [35], Mellet and Vasseur 
improve the results of Bresch,Desjardins by dealing with the case of a general pressure 
P{p) = ap^ with 7 > 1. 

In the case /t > 0, we remark then that the density belongs to L°°(0, 00, H^{Q,^)). Hence, 
in contrast to the non capillary case one can easily pass to the limit in the pressure term. 
However let us emphasize at this point that the above a priori bounds do not provide 
any L°° control on the density from below or from above. Indeed, even in dimension 
N = 2, functions are not necessarily locally bounded. Thus, vacuum patches are 
likely to form in the fluid in spite of the presence of capillary forces, which are expected 
to smooth out the density. It explains why it is so difficult to obtaining the existence 
of global strong solution in dimension N = 2. Indeed in order to prove the existence of 
global weak solution the main difficulty consists in dealing with the quadratic terms in 
gradient of the density appearing in the capillary tensor. Recently D. Bresch, B. Des- 
jardins and C-K. Lin in [8] got some global weak solutions for the isotherm Korteweg 
model with some specific viscosity coefficients. Indeed, they assume that p{p) = Cp 
with C > and \{p) = 0. By choosing these specific coefficients they obtain a gain of 
derivatives on the density p where p belongs to L'^{H^). It is easy at that time with this 
kind of estimate on the density p to get strong convergence on the term of capillarity. 
However a new difficulty appears concerning the loss of information on the gradient of 
u when vacuum existing (indeed the viscosity coefficients are degenerated). It becomes 
involved to pass to the limit in the term pnUn <^ Un- That's why the solutions of D. 
Bresch, B. Desjardins and C-K. Lin require some specific test functions which depend on 
the density p. 

In [26| . we improve the results of R. Danchin and B. Desjardins in [13] and D. Bresch, 
B. Desjardins an C-K. Lin in [8] by showing the existence of global weak solution with 
small initial data in the energy space for specific choices of the capillary coefficients and 
with general viscosity coefficient. Comparing with the results of |8], we get global weak 
solutions with general test function cp £ C^(r2^)not depending on the density p. In 
fact we have extracted of the structure of capillarity term a new energy inequality using 
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fractionary derivative which allows a gain of derivative on the density p. 
1.1 Derivation of the models 

We are going to prove that we can derive new entropy estimates when we choose specific 
coefficients for the viscosity and the capillarity. In the sequel we will consider the following 
physical coefficients: 

pip) = pp and k(p) = — , 

P 

with p,K > 0. By computation, we show that in this case, we obtain (see the appendix 
for more details): 

divi^ = Kdiv{pVVlnp) = Kdiv {pV{D In p)). 

It is now clear that we have to work with the new variable v = u + -V In p. We now want 
to rewrite system (jl.ip in terms of the variables (In p, v) . We have then by considering 
the mass equation: 

pdt In p + pu ■ VV lnp + pV In p ■ Vn* + pVdwu = 0. 

Then we obtain the following new system when a = ^: 

dtp + div{pv) - -Ap = 0, 

P (1.4) 

pdfV + pu ■ Vv — div{ppVv) + VP(p) = 0, 

When a = and k = p^, we obtain the following simplified model: 

dtp + d\v{pv) - -Ap = 0, 

P (1.5) 
pdtv + pu ■ Vv — div{ppVv) + VP(p) = 0, 

For more details on the computation, we refer to the appendix. Our goal is now to prove 
new entropy inequalities for these two systems and to obtain the existence of global weak 
solution for these two systems. We also prove the existence of global weak solution for 
the shallow-water system, indeed A. Mellet and A. Vasseur in [33j have obtained the 
stability of global weak solutions for shallow-water system but to get the existence of 
global weak solution, it remains to construct approximate solutions of the shallow water 
system conserving the entropy inequalities used in [33| . This is exactly what we will show 
by studying the system ()1.4p . 

2 Notations and main result 

2.1 Existence of global weak solution for Korteweg system 

We say that {p, v) is a weak solution of ()1.4p on [0, T] x i}, which the following initial 
conditions 

P/t=o = Po > 0, pu/t=o = mo. (2.6) 



4 



with: 

PQ £ L^{n) n L\n), ^vinpo e L\n), po > o, 

^ (2.7) 
^/PqVq £ L^(r2), Pq'^^ vq £ L^^'^(ri) for somme small 6. 



if 

• p£ L^{L^n) n LT(r?)), ^ G L^(L2(f))), 

• ^Vt; G L2((0,r) X n), 

with /) > and {p,^v) satisfying in V' {[0,T] x J]^): 

dtp + div(^^u) - -Ap = 0, 
p(0,x) = po{x). 

and if the following equality holds for all (p{t, x) smooth test function with compact 
support such that ^p(T, •) = 0: 

/ {pv)o ■ f{0,-)dx + / / ^{y/pv)dt^p + ® ^v: V ipdx 
Jn Jo JQ ^2 8) 

+ / / p'^divc^- < ppVv,Vip >= 0, 
Jo Jn 

where we give sense to the diffusion terms by rewriting him according to ^/p and ^/pv: 

< pVv,Vf >= — J y^{^Vj)dii(pjdx dt — J 2y/pvjdiy/pdi(pjdx dt 
Similarly we have exactly the same type of definition for system p.5|) . 

< pVvjVif >= — J ^/p{^/pvj)^iHPjdxdt — j 2^/pvj^i^/p^npjdx dt. 

Main results 

We obtain in this paper the existence of global weak solutions (more exactly the stability 
of global weak solutions) for systems ()1.4p and (jl.Sp . For system ()1.4p we obtain the 
following first theorem. 

Theorem 2.1 Assume that 7 > 1. Let {pn-,Vn) with Vn = Un + ^Vlog/3„ he a sequence 



of weak solutions of system \1.4\ satisfying entropy inequalities 113.18]) and I13.19\} . with 
initial data: 

{Pn)/t=0 = Poi^) (^^d {pnVn)/t=0 = PoVo{x) 

where Pq and Vq such that: 

Po>0, pI^Po inL\n), p^v^ ^ povo inL^Q), (2.9) 
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and satisfy the following bounds (with C constant independent on n): 

I {Po^^ + ^) <C, f ;^|VlnpSpdx < C, (2.10) 
J ft ^ 7 - J- Jn Pn 



and: 

,,n|2+a 
„n 1^0 



pV-^^dx < C, (2.11) 

Then, up to a subsequence, {pn^ \/PnVn^ \fPnUn) converges strongly to a weak solution 
(p, y^w,y^u) of |j.^[ ) satisfying entropy inequalities i3.18\) and 113. 19\} (the density pn 

3 

converges strongly in C^{{0,T), Lj^^^{Q)), ^Jf)^Vn converges strongly in L'^{0,T, Lf^^) and 
the momentum nin = Pn^n converges strongly in L^{0,T, Lj^^{Q)), for any T > 0). 

Remark 1 The proof will be strongly inspired of the works of A. Mellet and A. Vasseur 
in m 



We now obtain exactly thesame result for system (jl.Sp . 

Theorem 2.2 Assume that we have a sequence {pn,Vn) withvn = Un + -'V log pn of weak 
solutions of system U.5]) satisfying entropy inequalities \S.2U\) and (E Jifwith initial data: 

{Pn)/t=0 = Poi^) 0."^^ {pnVn)/t=0 = PoVo{x) 

where Pq and Vq such that: 

Po>0, Po^Po tnL\n), p^v^^povo inL\n), (2.12) 
and satisfy the following bounds (with C constant independent on n): 

[ {p^o^-^ + ^XC, [ ^\Vlnp-,fdx<C, (2.13) 
JQ ^ 1 - ^ Jn Po 

and: 

I n|2+<5 

Po-^dx < C, (2.14) 

Then, up to a subsequence, {pm ^/Pn'Vn, ^/PnUn) Converges strongly to a weak solution 
{p, ^v, ^u) of 111.5]} satisfying entropy inequalities 113. 20\) and 113.21]) (the density pn 

3 

converges strongly in C^{{0,T), L^^^{Q)), ^fp^Vn converges strongly in L'^{0,T, Lf^^) and 
the momentum rUn = PnVn converges strongly in L^{0,T, L}^^{n)), for any T > 0). 



I 

Jn 



In the specific case of the system ()1.5p we obtain new blow-up criterion in the case of 
the torus which improves the results in [T7]. In the sequel, we will set m = pv and 
q' = p — 1. We have then the following theorem: 

Theorem 2.3 Let P be a suitably smooth function of the density and 1 < p < +oo. Let 
tuq G Bp^oo with e > and q^ G Bp^oo such that po > c> 0. 
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There exists then a positive time T such that system U.5\) has a unique solution {q ,m) 
with p bounded away from and: 

~ —+e' ~ —+2+e' ~ —~l+e' ~ — + l+e' 

i^c{%nBi^ )nL}riB;^, ), mGC{[o,T];B;, ) n l\[o,t], b;^, ). 

We assume now that P{p) = ap with a > 0. Furthermore if vq G , pQ € Bp ,^ for any 
1 <p < +00, ^l{|po|<5} e ^°°(T^) n Li(T^) (with < S < I) and the initial data are 
in the energy space, it means: 

^/p^Uo £ L^, V^/p^ G and U{po) G L^. 

(we refer to the proposition \3. 1\ for the definition ofH) we can then extend the solution 
beyond (0, T) if: 

V £ LUL'^n^)) with 1 + — = i and l<p< +oo, (2.15) 

p 2q 2 

or if for any e > arbitrary small: 

^Mp\<S}^Lf{L\T'')). (2.16) 



Remark 2 This result is really to consider as a Prodi-Serrin theorem on the effective 
velocity v. In terms of blow-up condition, he improves widely J17j.In fact the second 
condition could be improved as follows (we refer to the proof of theorem \2. S\) : 

Remark 3 We could probably extend this previous result for more general pressure terms. 
However it would requires additional informations on the integrability of the density or 
of the vacuum, i.e ^l{|p|<<5}- 

It would be also possible to deal with the euclidian space (it does not change a lot, 
except that sometime, we shall use interpolation so that we would need of additional 
information in terms of low frequencies. 

3 New entropies 

3.0.1 Entropy for the system (|1.4p 

We now want to establish new entropy inequaUty for system ()1.4p and ()1.5p . We now 
obtain the following proposition when {p,u) are exact solutions of system (jl.4p . 
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Proposition 3.1 Assume that {p,u) are exact solutions of system (fi.^[ j with P{p) = ap^ 
fr > then for all t > 0; 

/ [p\u\'^{t,x) + K\Vy/p\'^{t,x) +U{p){t,x)]dx + [ [ V In p ■ V p^f dxdt 
Jn Jo Jn 

+ (m-7) / / p\yu\'^dxdt + j I I p\Du\'^dxdt + K / / p{dij\n pf{t,x)dtdx 



Jn Jo Jn Jo Jn 



<C{[ {po\vo\Hx) + U{poix)) + K\V^fit, x) dx). 
Jn« 

(3.18) 

with li{s) = sp j^^dz. 

Proof: We now want to obtain this new entropy by two different ways, one wliicli is very 
direct and an other using the BD entropy (see Q)- 

Direct energy inequality via system (II. 4p 

When we multiply the momentum equation in ()1.4p by we obtain: 

{p{t,x)\v{t,x)\'^ +I{{p){t,x))dx+ I I [pp{t,x)\Vv\'^{t,x) 
n JO Jn 

+ -P"ip)\Vpf{t,x))dtdx <C{[ {po{x)\vo{x)f + U{po){x))dx). 

Jn 

By the previous inequality and ()1.3p we obtain the desired result. 
BD entropy via system (jl.ip 

By an other way we obtain the new BD entropy of proposition ()3.ip by multiplying the 
momentum equation by V In p. Indeed the only difference with the classical BD entropy 
(see [Z]) is the extra capillarity terms and the term ^div (pVu*) in (jl.ip . We have only 
to treat these two extra terms. 

f [ dwK -V In pdxdt = K [ [ div{pVV In p) -V In pdxdt, 
Jo Jn Jo Jn 

= —K / / p{dij In p)'^dtdx. 
Jo Jn 

rt r ,^ rt 



/ / — div(pVn*) • V Inpdtdx = — / / —pdjUidij In pdtdx, 
Jo Jn ^ Jo Jn I'' 



K , 1 
< - — 

- p^2e 



[ [ p{dijlnpfdtdx + ^ f [ plVu'^fdx. 
Jo Jn 2 Jo Jn 



By bootstrap we conclude. For more details on the BD entropy in our case, we refer to 
the appendix. □ 

The following proposition comes from |34j. 
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Proposition 3.2 Smooth solutions of system satisfy the following inequality when 

P{p) = ap^ with 7 > 1/ 

(/j.--f)A..)'^(X*P..)i. (3.1.) 

/or 5 G (0,2). 

Proof: The proof follows exactly the same lines than in [33j . 
3.0.2 Entropy for the system (jl.Sp 

By proceeding similarly we obtain the following propositions for system (jl.Sp . 



Proposition 3.3 Assume that (p, u) are exact solutions of system ( ti.5|) then for all 
t > 0; 

/ {p\v\'^{t,x) + Ii{p{t,x))) dx + I I V In p ■ V p"' dxdt 

Jq Jo Jq ^2 20) 

+ K p{dij In pf {t,x)dtdx <C{ {po\vof{x)+U{po{x)))dx). 

Jo Jn Jn 

Proof: The proof follows the same lines as in proposition 13.11 □ 
For the following proposition, we refer to [33]. 

Proposition 3.4 The smooth solutions of system \1.5\ satisfy the following inequality 
when P{p) = ap^ with 7 > l: 

2 

d I pH^^^ I p\^\&\yy\2dx< ( [ (/,27-l-f)2^dxy''( / plvl^dx)'^, (3.21) 

for 6 e (0,2). 

4 Littlewood-Paley theory and Besov spaces 

Throughout the paper, C stands for a constant whose exact meaning depends on the 
context. The notation A < B means that A < CB. For all Banach space X, we 
denote by C([0,T],X) the set of continuous functions on [0,T] with values in X. For 
p G [1, +00], the notation LP{Q, T, X) or L^(X) stands for the set of measurable functions 
on (0,T) with values in X such that t — )• ||/(i)llx belongs to LP{0,T). Littlewood-Paley 
decomposition corresponds to a dyadic decomposition of the space in Fourier variables. 
Let if G C°°(M^), supported in the sheh C = e ft^/l < \^\ < §} and x G C~(M^) 
supported in the ball 5(0, |). (p and x are valued in [0,1]. We set = (0,27r)^ and 
= (Z/1)^ the dual lattice associated to T^. We decompose now u G S'{T^) into 
Fourier series: 



ui^) = 2^ u^e'P'^ with 2^ = -~ J ^e "^■yu{y)dy. 
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Denoting; 



one can now define the periodic dyadic blocks as: 



and the low frequency cutt-off: 



p<q-\ 



It is obvious that: 



feez 

This decomposition is called non-homogeneous Littlewood-Palcy decomposition. 
Furthermore we have the following proposition where C = i?(0, |) + C 

Proposition 4.5 



\k-]^\>2 =^ supp(^(2~''-) n supp¥:'(2-^ •) = 0, (4.22) 
k>l =^ suppx n supp<^(2-'=-) = 0, (4.23) 
I A; - A;' I > 5 =^ 2^'cf\2^C = $. (4.24) 

4.1 Non homogeneous Besov spaces and first properties 
Definition 4.1 For s G fi, p G [1, +oo], q G [1, +oo], anc? u G (S'(f2-^) we set.- 

ll^lls|.. = (E(2'llAn||L.)^)'. 

T/ie non homogeneous Besov space Bp ^ is the set of temperate distribution u such that 

||n|| OS < +00. 

p,<i 

Remark 4 The above definition is a natural generalization of the nonhomogeneous Sobolev 
and Rdlder spaces: one can show that B^ ^ is the nonhomogeneous Rdlder space and 
that B2 2 is the nonhomogeneous space . 

Proposition 4.6 The following properties holds: 

1. there exists a constant universal C such that: 

C^^WuWb' < ||V?i||n3-l < CIImII Rs . 
" " p,r — II ll-t>p,r — II II p,r 

2. Ifpi < p2 and n < r2 then B^^^,^ ^ B^'f^^^'^'-^'^'K 
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3. -Bp^n ^ Bp J. if s' > s or if s = s and ri < r. 

Let now recall a few product laws in Besov spaces coming directly from the paradifferen- 
tial calculus of J-M. Bony (see [1]) and rewrite on a generalized form in [1] by H. Abidi 
and M. Paicu (in this article the results are written in the case of homogeneous sapces 
but it can easily generalize for the nonhomogeneous Besov spaces). 

Proposition 4.7 We have the following laws of product: 

• For all s G {p,r) G [l,+oo]^ we have: 

Ww^Wb^p^, < C{\\u\\loo\\v\\b=^,. + \\v\\l^\\u\\b=^^^) . (4.25) 

• Let (p,Pi,P2,r, Ai, A2) G [l,+oo]2 such that:^ < + j^, Pi < X2, P2 < Ai, ^ < 
^ + and - < 7^ + 1^. We have then the following inequalities: 

ifsi + S2 + Nm/{0, 1 -± - ±) > 0, si + < ^ and S2 + < then: 

\\uv\\ < \W\\b;iJv\\sS2^^, (4.26) 

when si + ^^=^ (resp s2 + ^ = ^)we replace \H\ b;iJv\\ b;i^ (^^^^P Mb^iJ 

by \\u\\B;iJ\v\\B;ir (''^'P MB;i^nL^)' + £ = f "^'^ ^2 + f = f w^e take 
r = 1. 

Ifs, + S2 = 0, «ie(f -f] and + < 1 then: 
If \s\ < — for p > 2 and —^■ < s < — else, we have: 

J \ \ p J f — p p ! 

WuvWb" < C'llnllRs n . (4.28) 

II ll-^p,r — II ll-^p,rll oT" r 

Remark 5 In the sequel p will be either pi or p2 and in this case j — ^ ~ ^ 'if Pi ^ P2, 
''^'Pj = V2-jl ifP2<Pi- 

Corollary 1 Let r G [1, +oo], 1 <p <pi < +oo and s such that: 

. s G (-^,^) + ^ < 1, 
^ PI pi ' p pi — 

, s G (-^ + iV(i + ^ - 1), ^) i/ i + ^ > 1, 

\ pi ^p pi ' ^ pi ' ■' p pi ' 

N_ 

then we have if u £ Bp ,, and v G Bpl^oo n -L°° : 

\\uv\\b' < C||u||rs \\v\\ n 

II II -^p.r — II II -^p.r II II ■:pr- 

^p-l,oo' 1-^ 
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The study of non stationary PDE's requires space of type LP{0,T, X) for appropriate 
Banach spaces X. In our case, we expect X to be a Besov space, so that it is natural to 
locahze the equation through Littlewood-Payley decomposition. But, in doing so, we ob- 
tain bounds in spaces which are not type LP{0, T, X) (except if r = p). We are now going 
to define the spaces of Chemin-Lerner in which we will work, which are a refinement of the 
spaces L^(i?p ,.). 



Definition 4.2 Let p € [l,+oo], T G [l,+oo] and si G Q. We set: 

ll^lli;^(B-) = (E2''"l|A/n(t)||l.(^.))^ 



We then define the space Lj,(S*i^) as the set of temperate distribution u over (0, T) x Q,^ 
such that \\u\\rp ,ns-^ , < +00. 

We set CriBp^r) = ^ri^p^r) <^ C{[0,T], Bp^^.). Let us emphasize that, according to 
Minkowski inequality, we have: 



Remark 6 It is easy to generalize proposition \4- 7\ to Lj,{Bp^^) spaces. The indices si, 
p, r behave just as in the stationary case whereas the time exponent p behaves according 
to Holder inequality. 

In the sequel we will need of composition lemma in Lj,{Bp j.) spaces. 
Lemma 1 Let s > 0, {p,r) £ [1, +00] and u £ L^{Bp ,.) D Lf{L°°). 

1. Let F G wl^l'^'^'°°{n^) such that F(0) = 0. Then F{u) G L^{B^,.). More precisely 
there exists a function C depending only on s, p, r, N and F such that: 

II^MIIZ^(B|,,,) ^ ^(ll^lli§?{i°°))ll^llL^(B|,J- 

2. Let F G Wl'}^^'^{il^) such that F{0) = 0. Then F{u) - F' {0)u G L^{B^^r)- More 
precisely there exists a function C depending only on s, p, r, N and F such that: 

\\F{u) - F'(0)tx||jp < C7(h|Uoo(i^))||n|||, 
Let us now give some estimates for the heat equation: 

Proposition 4.8 Let s G il, {p-,r) G [1, +00]^ and 1 < /O2 ^ Pi ^ +00. Assume that 
uq G Bp j. and f G L'^{Bp^r^^'^^'^^). Let u be a solution of: 

dtu — pAu = f 
Ut=0 = Uo . 

Then there exists C > depending only on N, p, pi and p2 such that: 

^ c(II^oI|b|,, + /^^"' 11/11 jP2(^.-2+2/p,^) . 

If in addition r is finite then u belongs to C([0, T], ,.). 
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5 Proof of the theorems 12.11 and 12.2 



We now present the proof of theorem 12.11 To begin with, we need to make precise the 
assumptions on the initial data. 

Initial data: 

We recall that the initial data must satisfy ()2.10p . and ()2.1ip to make use of all the 
inequalities presented in the previous section: 

• is bounded in L^{n) n L'^(O), p^>0 a.e in Q, 

• /OokoP is bounded in L^{^), 

• 111 Po is bounded in ($7), 
Po\vl^\^^^ is bounded in L\n). 



With those assumptions, and using the entropies ()3.18p . (|3.19p and the mass equation, 
we deduce the following energy estimates, which we shall give uniform bounds: 



ll\/Pnll-L°°{(0,T),L2(Q)) < C, 
ll/'n||L°o{(0,T),LT{n)) < C, 
\WPndijlnpn\\L2(^{0,T)xn) < C, 
\\\^^Un\\L2{(o,T)xn < C, 

and 

\\Pn V/9„||i2((o,T)xn) < C, 
\\Pn\Vn\ \\L°°{{0,T),L^n)) < C. 

In view of previous inequalities, the bounds ()5.29p and (|5.30p yields: 

\\^/p^Vn\\L2{(0,T)xn) < C, 

l|V\/Pn|lL°°(0,T;L2(n) < C, 
X 

ll^Pn llL2((o,T)xn) < C, 
\\Pn\Vn\'^~^^\\L'^({0,T),L^(n)) < C. 



(5.29) 



(5.30) 



(5.31) 



The proof of theorem 12.11 will be derived in 5 steps and follows the proof of [33j . In the 
first two steps, we show the convergence of the density and the pressure (note that the 
convergence of the density is straightforward here). The key argument of the proof is 
presented in the third step: as y/p^Vn is bounded in a space better than L°°(0,T; L^(f])), 
it will allow to give the convergence of the momentum (step 4) and finally the strong 
convergence of y/p^Vn ^/IhiUn in L^((0, T) x 0) (step 5). In fact by considering the 
effective velocity Vn we do not need extra information on the density pn to deal with the 
capillarity terms (as it is the case in [3126]). Furthermore at the difference with [8] we are 
able to treat the momentum term without using test functions depending of the density 
(or in an other way to add a cold pressure). The last step address the convergence of the 
diffusion terms and of ^fp^Vn to y^n + ^V/>. 
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Step 1: Convergence of ^fp^ 

Lemma 2 We have the following properties: 

1. is bounded in L°°{0,T; H^{n)) 

2. dt^/p;^ is bounded in L'^{<d,T;H-^{Vi)). 

As a consequence up to a subsequence, converge almost everywhere and strongly in 
L'^iO,T;Ll^{n)). We write: 

-y/p^ — )• ^/p a.e and Lfg^{{0,T) x strong. 

3 

Moreover, pn converges to p in (7*^(0, T; L^^^^(il)). 

Proof: The second estimate in (j5.3ip . together with the conservation of mass HpnllLi = 
||/9n,o||Li gives the L'^{0,T; H^{^})) bound. Next as we have: 

dt\^p^ = -^VPndivu„ - Un ■ V^/Jhi 

thanks to Aubin's Lions Lemma, we conclude to the strong convergence in L^g^((0, T) x Vt) 
of Up to a subsequence we obtain that pn converges to p a.e. 

Sobolev embedding implies that is bounded in L°°(0, T; L'?(r2)) for q G [2,+oo[ 
if = 2 and q G [2, 6] if = 3. In either cases we deduce that pn is bounded in 
L°°(0,T;L3(17)), and therefore: 

3 

pnUn = y/Pn^/PnUn is bounded in L°°(0, T; L2 (ri)). 

The continuity equation thus yields dtPn bounded in L°°(0, T; VF~^'2 (17)). Moreover 
since Vp„ = 2y^Vy^, we also have that Vp„, is bounded in L°°(0, T; La (J7)), hence 

3 

the compactness of pn in C{[Q,T]\Lf^^{Q)). 
Step 2: Convergence of the pressure 

Our goal now is to prove that the pressure term is enough integrable in L^+'^((0, T) x 17). 
Indeed it will be sufficient to prove that P{pn) converges to P{p) in L\^^{^, T) x 17) since 
we know from the previous lemma that pn converges a.e to p. 

Lemma 3 The pressure pZ is bounded in Ls ((0, T) x 17) when N = 3 and L''((0, T) x 17) 
for all r G [1, 2[ when N = 2. In particular, pn converges to p^ strongly in L[^^((0, T) x 17). 

Proof: Inequalities ffOT]) and (OH yield pi G L'^{0,T; H'^). 

When Af = 2, we deduce p| G L2(0, T; L'?(17^)) for all q G [2,oo[. So pZ is bounded 
in L^(0,T] LP) n L°°{L^) for all p G [l,+oo[, hence by interpolation pZ is bounded in 
U{{Q,T) X 17) for all r G [1,2[. 
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When = 3, we get by Sobolev embedding that pn is bounded in V-{{),T;L'^{^)). As 
Pn is also in L°°((0, T), L-^(r2)), by interpolation we have: 

2 3 
IITII <^II7I|5 IITII'5 

"^"llif ({0,T)xQ) - ll^«llL°°{0,T;Ll(n))ll^™llLl{0,T;L3(n)) 

Hence pn is bounded in L3((0,T) x 0). 

Since we already now that pZ, converges almost everywhere to p^ , we can conclude to the 
strong convergence of pZ, in Lj^^{{0,T) x fl). 

Step 3: Convergence of the momentum 

Lemma 4 Up to a subsequence, the momentum rUn = PnVn converges strongly in L^(0, T; L-f^^{Q)) 
to some m{x,t) for all p G [1, |). In particular: 

PnVn rn almost everywhere {x,t) G (0,T) x Q. 

Note that we can already define v{t,x) = "^j^l'^^^ outside the vacuum set {p{t,x) = 0}, but 
we do not know yet whether m{t, x) is zero on the vacuum set. 

Proof: We have: 

Pn^n — Pn\l Pn^ni 

where ^ is bounded in L°°(0, T; for q G [2,+oo[ if iV = 2 and ^ G [2,6] if 

= 3. Since ^ffhivn is bounded in L°^(0, T; we deduce that: 

PnVn is bounded in L°°(0, T; L«(17)) for ah gG[l,^]. 

Next we have: 

diiPnvi) = PndiVnj + u{diPn 

= \/%^/p^.divl + ^^j^vndi^fp:;^. 

Using (j5.3ip the second term is bounded in L°°(0,T; while the first term is 

bounded in ^^(o, T; for all g G [1, |]. Then: 

^{PnVn) is bounded in L^(0, T; 

In particular we have: 

PnVn is bounded in ^^(0, T; W'^'^(O)). 

It remains to show that for every compact set X G fi, we have: 

dt{pnVn) is bounded in L^{{),T;W~'^^^K)). (5.32) 

Indeed ()5.32p and Aubin's Lemma (see Lions [31j), yields the compactness of pnVn in 
L\0,T;LP{K)) for ah p G [1, |). 
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To prove ()5.32p . we use the momentum equation of system ()1.4p . first we observe that 
using lemma E] and (|5.31|) : 

VpleL^{Q,T-W-^^\K)). 

So we only have to check that the terms V{pnVvn), V{PnVvn) are bounded in L°°{0, T; W~ 
To that purpose, we write: 

/OnVf„ = S/{pnVn) - Vn'Vpn, (5.33) 

(and we proceed similarly with the other two terms). The second term in ()5.33p reads: 

which is bounded in L°°(0, T; L'^{n)) thanks to (iOOD and (iOT]) . The first term in (iOHIl 
can be rewritten: 

V[/0„V„] = V[^/p:^{^/p:^Vn)], 

which is bounded in L°°(0,T; W~^'^{Q,)) thanks to lemma|21 
We deduce that pnD{vn) is bounded in: 

L°°(0,r;VF-^'t(K) + Li(K)), 

and since L^{K) C VF"^'t(/C) and W~'^^i{K) C W'^^^{K) we conclude that PnD{vn) 
and /9„ are bounded in L°^(0, T; 's (X)), which conclude the proof of lemma [H 

Step 4: Convergence of ^fp^Un ® y/Pn^n 

Lemma 5 The quantity ^j^Vn converges strongly in Lf^^{{0,T) x Q,) to ^ (defined to 
be zero when m = 0). 

In particular, we have m{t,x) = a.e on {p{t,x) = 0} and there exists a function v{t,x) 
such that m{t,x) = p{t,x)v{t,x) and: 

3 

PnVn pv strongly in L^(0, T; Lf^^(ri)), p G [1, -), 

y/p^Vn ^/pv strongly in Lf^^{{0,T) x Q), 
(note that v is not uniquely defined on the vacuum set {p{t,x) = 0}). 
Proof: First of all, since is bounded in L°°(0, T; L2(17)), Fatou's lemma yields: 

liminf — —dx < +oo. 

Pn 

In particular, we have m{t, x) = a.e. in {p{t, x) = 0}. So if we define the limit velocity 
by v{t,x) by setting v{t,x) = when p{t,x) ^ and v(t,x) = when p{t,x) = 0, 

we have: 

m{t,x) = p{t,x)v{t,x) 
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and 



— dx = J p\v\'^dx < +00. 



Moreover, Fatou's lemma yields: 

J p\v\'^^^dx < J lira inf pn\vn\'^^^dx 

< liminf / pn\vn\'^^^dx, 



and so p\v\'^+^ is in L°°{0,T; L^{n)). 

Next, since rUn and pn converge almost everywhere, it is readily seen that in {pit, x) 7^ 0}, 
yJ'p^Vn = ^= converges almost everywhere to y/pu = Moreover, we have: 

\fP^'"n'^{\vn\<M} \fpv'^{\v\<M} almost everywhere. (5.34) 

As a matter of fact, the convergence holds almost everywhere in {p{t,x) 7^ 0}, and in 
{p{t,x) = 0}, we have y/p^Vnl{\v^\<M) < M^/p^^ 0. 

We are now in position to complete the proof of lemma O For M > we cut the 
norm as follows: 



PnVn-\/pv\ dxdt < / \y/p^Vnl{\y„\<M} - Vpvl{\v\<M}\ dxdt 



+ 2y \y/p^Vnl{\v„\>M}\'^dxdt + 2 J \^vl{i^i>M}\'^dxdt, 

It is obvious that y^Wnl{|j;„|<M} is bounded uniformly in L°°(0, T; L^{^)) and by ()5.34p . 
we obtain that y/p^Vnl{\y„\<M} converges strongly to ^/pvl{\v\\<M} in Lf„^((0,r) x 0): 



^Vni{\v„\<M} - ^/pv'^{\v\<M}\'^dxdt 0. 



(5.35) 



Finally, we write: 



and 



^Vnl{\v„\>M}\'^dxdt < 



p\Vr 



\2+8 



dxdt. 



(5.36) 
(5.37) 



J \VPvl{i^i>M}\'^dxdt < ^ j p\v\^^^dxdt, 
Putting together ()5.35p . ()5.36p and (|5.37p . we deduce: 

/C 
\^/p^Un - y/pu\'^dxdt < —J-, 
IVl 

for all M > 0, and the lemma follows by taking M — )• +00. 

Lemma 6 ^/p^Un fX" yfPnVn converges to the distribution sense to yfp^Un ® y/PnVn- 

Proof: It suffices to observe that ^fp^Un converges weakly in ^^(0, T, ^^((7)) to 
and y/p^Vn converges strongly in Lf^^{{0,T) x 17) to y^f. We have then for all (/? G 

c^{{o,T)xny. 



Jn 



{y/P^Un (8) y/p^.Vn)(pdxdt 



JQ 



{y/pu (g) y/pv)(pdxdt. 
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Step 5: Convergence of the diffusion terms 

Lemma 7 We have: ^ 

Pn^Vn — )• pVv in V , 

Pn^Vn />*Vf in V' , 
Proof: Let (/> be a test function, then: 

Pn^Vn4>dxdt = — PnVn^ (pdxdt + / VnV Pn 4>dxdt 



^i{y/p^Vn)V 4>dxdt + j ^/pl^VnV {2./p^)(j)dxdt. 



Thanks to lemma [21 we know that ^/fhi is bounded in L°°(0, T; L^q^(O)) and that ^yp^ 
converges almost everywhere to ^/p (defined to be zero on the vacuum set). Therefore it 
converges strongly in Lf^^{(0,T) x il). This point is enough to prove the convergence of 
the first term. 

Next as converges strongly to ^ in Lf^^{{0,T) x $7), it follows that: 

^ V^P ^L((0, T) X 0) - weak. 
And we conclude for the second term as ^/p^Vn converges strongly in L^^^((0, T) x Q,). 

Convergence of ^/p^Vn to + ^Vy^ 
Lemma 8 We have the following property. 

^ r- ^ 

Proof: As: 

Ap„ = v7^(\/p^Alnp„) + 4|Vv^'^ 



n| ) 

we can show by (fOgj) that Ap„ is uniformly bounded in Lf{L^{n)) + L^L2{n)). By 
Aubin Lions theorem, we deduce that Vp„ converges a.e to V/5. It concludes the proof 
of the lemma. 

5.1 Stability for theorem [2721 

The proof follows the same line than for theorem 12.11 and than [3^3j . 



6 Proof of theorem 12.3 



For the result of existence of strong solution in the theorem 12.31 we refer to [T7| . We 
now want to concentrate on the result of blow-up of the theorem ??. To do this, we 
shall point out new regularizing effects and gain of integrability on the density and the 
effective velocity v. Indeed the key of the proof is to obtain enough integrability on v to 
derive regularizing effects on the density via the first equation of system 11.51 (which is a 
heat equation with a remainder term which depends on v, i.e div{pv)). 
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Gain of integrability on v 

As in |24J, we now want to obtain additional information on the integrability of and 

- 

more precisely we would like to show that p^v is in any L^f (LP(T^) with l<p< +00. 
To do it, we multiply the momentum equation of (jl.Sp by tjl^;!^"^ and integrate over 
we obtain then: 



- / pdt{\v\'P)dx + f pn-V(^)dx+ / p\v\P-^\Vv\^dx 

P JjN JjN p JjN 

+ (p — 2) / py^ VjVkdiVjdiVk\v\^ '^dx+ / \v\^ '^v-\7apdx = 0. 



+ 



-t 




By integration by parts we have: 



a\v\P-^vVap(t,x)dtdx = - [ [ diY{\v\P~'^v)ap{t,x)dtdx. 



(6.39) 



Next we observe that: 

^VjVkdiVjdiVk = ^(^VjdiVjf = ^ ^9i(|?;p). 

i,j,k i j * 

We get then as div{pu) = —dtp and by using (j6.38p : 

- / dtip\v\P)dx+ [ p\v\'P-^\Vv\'^dx + ip-2) [ p(y ^^{\v\'^f\v\P-^dx 

+ / \v\P~'^v ■ Vapdx = 0. 
We have then by integrating over (0, t) with < t < T: 

- / {p\v\P){t,x)dx+ [ [ p\vf~^\Vv\'^{t,x)dtdx 

P JjN Jq JjN 

{p-2) [' [ p{y2d^{\v\^f\vrHt,x)dtdx<- [ {po\vo\n{x)dx (6.40) 

Jo Jjf^ , P JjN 

\v\P^'^v ■ \7ap{t, x)dtdx\. 

We now want to take the sup of the previous estimate on (0,T), we have then: 

sup (- / {p\v\P){t,x)dx + [ [ p\v\P'^\Vv\^{s,x)dsdx 
te(o,T) P JT'v Jo Jt^ 

+ {p-2)f [ p{y^d,{\v\''f\v\'"\s,x)dtdx) <- I {po\vo\n{x)dx (6.41) 

Jo JjN i P JjN 

+ sup I / / \v\P~'^v ■ 'Vap{s,x)dtdx\. 
te(o,T) Jo JT^ 
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We have then: 

div{\v\P-^v) = \v\P-^div{v) + {p - 2)\v\p-\ ■ {v ■ Vv). 

We now want to bound sup^g^Q | Jq J^j^ \'^\^ ' ^CLp(t,x)cltdx\ in (|6.40p . To do this, 
we would Hke to apply Holder's inequalities, indeed we have: 

\p\v\^~ V ■ {v • Vv)\ < C\^/p\v\2''-Vv\ x|/)2 p|ti|2^-^|x|/)p|. 
By Holders inequalities we can prove that IpIvIp-'^v ■ {v ■ Vu)| is in L]-{L}{T^) with: 



\p\v\P ■ {v ■Vv)\dxds < sup (v^||\/p|'u| ^ '^^v\\i^2(^i^2(jn-^-^ 



sup 

t6(0,T) Jo Jf^ te{0,T) 



1 



< sup {Vi\\^\v\^'^Vv\\L2(^L2nN)) sup ( / p\v\\P{x)dx)'^'p\\pp\\Lc^^Lp^jN))). 
te(o,T) se(o,i) Jt'^ 

(6.42) 

Now we set: 

A{t) = (-[ {p\v\P){t,x)dx+ [ [ p\v\P'^\VvWs,x)dsdx 

P JjN Jq JjN 

+ ip-2) r / p{yd^(\vmvr\s,x)dtdx). 

Jo Jj^ ^ 

We have then from (|aIT]l and (j02|) : 

sup A{t)<C{l + VT{ sup A{t))^-p\\p\\l^ ) 

ie{0,T) t6(0,T) T >. V JJ 

1 

We deduce then that sup^g^o T) finite. In particular, we obtain that ppv is in any 

L^{LP{T^)) for 1< p < +00. 

Blow-up criterion 

We have then obtain that for any 1 < p < +oo, ppv belongs to L°°[LP). We now want 
to take into account this information to obtain regularizing effects on the density via the 
first equation of (jl.Sp (which is a heat equation): 

dtp - - A/9 = -d\\{pv). (6.43) 

Our goal is to transfer the information on the integrability of v (which is a subscaling 
estimate) on the density p. More precisely we have by proposition 14.81 for any 1 < p < 
+oo: 

Ml^iBl^) ^ C{\\po\\i^^B^^^^ + WpvWl^^Bo J. (6.44) 
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We now need to prove that pv is in L'^{Bp ,^). As ppv is in L'^{U'), we have then: 

\\{P '•)P''v\\l^{Lv{JN))<\\P''v\\l^(^Lv)\\p HlL5f(L-(T^))> 

< C||p''^llL-(LP)ll/5|li/(ioo(TriV))- 



By injecting (j6.45p in (j6.44p and using the fact that LP(T^) is embedded in B^,^, we 
obtain that: 



\\Ph^(Bl^) ^ ^^\\P^h^{Bl^) + \\P~'v\\l'^{Lv(T^))\\p\\l^\l^(jn^))- (6-46) 

By Besov embedding we know that Bp ^^) is embedded in L°°(T^) for p large enough. 
We obtain by (lOTl) that: 

Mz^iBl^) ^ C'(||po|lz^(Bi^^) + IIP''^llL-(LP(TiV))||p|li/(Bi_^))). (6.47) 
We conclude then by bootstrap that p is in L^{Bp ,^) for p large enough (in fact p > N). 

Estimates on ^ if we control j;r in Lf{L'^{T'^)) 

We now would like to obtain estimates on ^, to do this we will consider the first equation 
of p.5|) that we multiply by — ^ with p > 2: 

' ^*(7^)-:7;;^^(7^) + :7^^|v(4T)l^ = -idiv(H, (6-48) 



p-l '^pp~^' p{p-l) pP-^ p{p-l)^ p2~h" pP 
because we have: 



Our goal is now to integrate over (0,t) X the equality (fOHjl . we obtain then: 



p-l JjN pP '^{t,x) p{p-lY Jq Jtn p2-\{t,x) 

1 \ f ^ 

-j dx+ \ / / — <l\\{pv){t,x)dtdx\. 

(x) Jo Jtn pP 



1 

< 



(6.49) 



We now pass to the sup on (0, T), and we have: 

< I „_! dx+ sup (I / / ^div(/)v)(t,x)(it(ij;|). 

P-^ Jin \x) te(o,T) Jo Jtn pP 

By integration by parts we get: 

— d\v{pv){t,x)dtdx\ = \p'S^ I I — dipv^{t,x)dtdx\. 

jjv pP ' ^ Jo JjN pP 
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(6.50) 



We recall that: 

\\dipv%t,x) = \^dip\ X X \v%t,x). (6.51) 

In the sequel we only shall deal with the case = 3, the other cases are similar to treat. 
Now as on the left side of ()6.50p . we have information on i in LS?(L^(T^)) and on 

V( p"'_i ) in L^(L^(T^)) we can assume by Sobolev embedding and interpolation that: 

^g4(L'^(T^)) withi + A = i. (6.52) 

We have then: 

1 1 i-| 1 I 

_ 2 

< ,, 

p2^2 ~-l ^~ " p2^2 "-i ^" " p2^2 

(6.53) 

We now want to take advantage of the fact that belongs to L^{L^{T^)) for e arbitrary 
small. Indeed the idea is now to use interpolqtion results in order to absorb the term on 
the right side of ()6.50p by the left side. More precisely we have by ()6.52p that: 

e L^{L\T^)) and ^ G L^^'^^(L^('^)(T^)), (6.54) 



1 1—- 1 1 2 

"eTT 1 1 l/(L2 (T^v ) ) ( 1 1 ^ ( ^TT I L2 (L2 (T^ )) + ^1 1 ^ZT 1 1 L?? (L2 (T^^ )) ) " • 



with: 



_ II ^ iiIF^ 



Now by interpolation we can show that: 



(6.55) 



ll^llL^(L^(TiV)) < ll;^lli^%i(TiV))ll;^ll^(i^)^^,(^,^^^^^> (6-56) 

with i = ^ and i = (1 - 0) + ^ and < < 1. 
We check easily that: 

1 1 ^ 

I|-tIIl9,(l/9(t^)) = II £_i II ^2«. 2,fs ■ (6.57) 
P p2 2 l|^(lf^(t'V)) 

From ()6.55p . ()6.56p and ()6.57p . we obtain finally that: 

1 1 1 2.e 

II^ItII 2m, 2.13 ^ II77IIl2?(L1(T'^))" £-1 lli^(i''{T'V))' 

or: 



In . „ 1 „(i^T)(l-^) II 1 ||((i^T))'e 



< ||_IHp-i)'- ' II : ||Mp-i)' f6 591 



i"Lf*(LF^(TiV)) - V^lli??(iHT-))ll H„i llL^(L.(T-)) 
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We recall that: 

^ = * „d ^ = _--^ii- (6.60) 

p-l 9 p-1 {1 - 9){p - l)q + 2e9 ^ ' 

By using ()6.51|) . ()6.53p and ()6.59p we obtain by Holder'sinequelities: 



/ / -^div(pt;)(t,j;)(it(ix| < C||^i-Y5i/9||i2(i2(TriN 

Jo JjN pP p2 + 2 ^ 



P 

1 

^ II r II rkil + a-,) ,^ „n-L^„> /^^^r^^ llf^llrfel^rnl 

P 



(6.61) 



£-1 llLy'^+"i'(L'J(i+"2)(T^))ll'^llL^^(L'J^(T^))- 



with Qi > 0, 02 > arbitrary small, 9 depending on oi > and 02 (such that < 9 < 1 
and 9 goes to 1 when oi and 02 go to 0) and k^, k"^ such that: 

2 + A;(l + Qi) + ^ " ^ 2 + ^(1+^+^2-1- 

From ()6.53p . ()6.6ip and ()6.59p we have finaly: 

I / / —diY{pv){t,x)dtdx\<C{ sup B{t)Y^^^)^^^ 

Jo JT^ P^ t£(0,T) 



= (0,T) 

II II 

I pe IIl5?(L1(T'V))II'^IIl*^(L<j1(T'V))' 



(6.62) 



with: 



When (p — 1) > 2e we can absorb the term | Jq fj^ ^div(/9t')(t, x)dtdx\ because from the 
previous inequality, we have: 

( sup Bit)Y<C{l + { sup B{t)f^^T^^"')\\v\\^,^ ^ ), 

fG{0,T) te(0,T) T 1) 

with 1 + ( )^^ < 2. Indeed we can show easily that v is bounded in L^^ (L*?^ (T^)). 

To do this, we recall that p^v is bounded in L^(L^'(T^)) for p enough big. As we control 
in L^(Li(T^)) and that: 

1 1 
V = —ppv, 

pv 

V is then bounded in L5?(L2 (T^)) for p enough big and then v is bounded in (L*? (T^)). 
We have then proved that for any p > 0, ^ is bounded in L^(L^(T^)). 

Estimates on ^ if we control v in L^(L''(T^)) with | + ^ = ^ and 1 < p < +00 

In this case, it suffices only to bound | Jg jr^j^ ^(li'v{pv){t, x)dtdx\ as follows: 



^div(pw)(t, x)(it(ix| < C|| ^T^iPlL^ (L2(T^)) II II , ' ll^'llL^(L'/(T^))• 



23 



And by Sobolev embedding we have: 

I / / — diY{pv){t,x)dtdx\ < C( sup -B(i))^||f H^^p 
Jo Jt!^ te(o,T) ^ 



= (0,T) 

This term can be absorbed in (|6.5U|) because if T is small enough, then 11^11^^ (^gcriv)) 
is small enough to do a bootstrap. To obtain the result on a general interval (0, T) it 
suffices to decompose the interval in a sum of small intervals and to apply the previous 
idea by iteration. 

Control of Vlnp in L°°{B%^^^^) 

To do this it suffice to recall that is in Lf{L'^{T^)), and as ^ is in L^{L^{T^)) 

for all < a then by Sobolev embedding we can show easily that ^ is in L^{B^_^_^ ^) 

I fnr anv -n ^ 1 Yi'^r 

P,' 



Since we have seen that Vp is in L^{Bp ,^) for any p > 1, by paraproduct we obtain 



that: 



Vlnp e L°°(5° 



N+e,oo)- 



Furthermore we have obtained that ppv belongs to L°°{I/') for any 1 < p < +00 and 
that Vlnp G (^^+e,oo)- As u = V — -Vln/9, we obtain easily by following the same 
lines that u is in L^{B^^^ ^. To summarize what we have obtained, we have: 

ueL^iB%^,^^) and Inp G L??(i?i.+,,^). (6.63) 

Easily we can prove that if uq G -Bjv+e 00) ^^'^ Po ^ ^N+e 00 then the system ()1.5p has 
a strong solution (p, u) on (0, T ) with: 

C 

T > 



(l + ||uo||ijo +||ln/9o|lBi ; 



We refer to [T7] for the proof. This is an easy consequence of the fact that the initial 
data are choose subcritical. 

It means that there exists a time T > c > 0, where c depends only on the physical 
coefficients and of subcritical initial data. We can construct by theorem ?? a solution 
(pi, ui) on (T — a,T — a + T ) with initial data {p{T — a),u(T — a)) (here a < T ). The 
only difficulty is to prove that on (T — a, T) we have: 

(pi,ni) = {p,u). 

To do this, it suffices only to use the uniqueness part of theorem 1 in [17'. It concludes 
the proof of theorem 12.31 □ 

7 Appendix 

In this appendix, we only want to detail the computation on the Korteweg tensor. 
Lemma 9 

divK = ft;div(pVVlnp) = Kdiv (pZ) ( V In p)). 
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Proof: By calculus, we obtain then: 

I 

(divK)j = (VAp - div{-Vp Vp)) ., 

= d,Ap--Apd,p-^d,\Vp\^ + ^\Vp\^djp, 
p zp p 

Next we have: 

Ap = pAlnp+ -|Vpp. 



(7.64) 



We have then: 



1 ■ " «... |2n a i„ „a „ ^ lv7„|2f 



(7.65) 



djAp- -Apdjp = dj{pAlnp+-\Vp\'') - Alnpdjp- -^\Vp\^djp, 

= pd,A\np+U,{\Vp\^) - ^\Vp\^d,p. 

Putting the expression of (|7.65|) in (|7.64p . we obtain: 

(divi^), = djAp+^dji\Vp\') - ^\Vp\'d,p. (7.66) 
zp p 



Next by calculus, we have: 



^a,(|Vp|2) - ^\Vp\'^jp = Y,i^^lnp^^,p-i^^lnpf^jp), 
P 

= diln p pdi J In p, (7.67) 



= |V(|lnp|2),-. 

Finally by using ()7.67p and ()7.66p . we obtain: 

divK = p{VA{lnp) + |v(|V Inpp)). 

We now want to prove that we can rewrite ([7]) under the form of a viscosity tensor. To 
see this, we have: 

div(pV(Vln/>))j = ^diipdijlnp), 

i 

= X] \-^iP^-i-i P + P^^iij P] ' 
i 

= p{AV In p)j + ^ pdi In pdjdi In p), 

i 

= p(AVlnp), + |(V(|Vlnp|2)),., 
= divi^. 

We have then: 

diviv: = ft;div(pVVlnp) = Kdiv(pi:»(Vlnp)). 
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